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Description 

Field of the invention 

5 [0001] The invention generally relates to techniques for analyzing biological samples, such as DNA or RNA samples, 
and in particular to techniques for analyzing the output of a hyt>ridized biochip to which the sample has been applied. 

Background of the Invention 

10 [0002] A variety of techniques have been developed to analyze DNA, RNA samples or other biological samples to 
identify diseases, mutations, or other conditions present within a patient providing the sample. Such techniques may 
determine, for example, whether the patient has any particular disease such as cancer or AIDS, or has a genetic 
predisposition toward the disease. 

[0003] One particularty promising technique for analyzing biological samples uses a DNA biochip (or microarray) 

IS which generates a hybridization pattem representative of the characteristics of the DNA within the sample. Briefly, a 
DNA microarray includes a rectangular array of single stranded DNA fragments. Each element within the array includes 
millions of copies of identical single stranded strips of DNA containing specific sequences of bases. A different fragment 
of DNA may be provided at each different element of the array. In other words, location (1,1) contains a different single 
stranded fragment of DNA than location (1,2) which also differs from location (1,3) etc. 

20 [0004] A DNA sample to be analyzed is first fragmented into Individual single stranded sequences with each se- 
quence being tagged with a fluorescent molecule. The fragments are applied to the microarray where each fragment 
bonds only with matching DNA fragments already embedded on the microarray. Fragments which do not match any 
of the elements of the microarray simply do not bond at any of the sites of the microarray and are discarded. Thus, 
only those microarray locations containing fragments that match fragments within the DNA sample will receive the 

25 fluorescent molecules. Typically, a fluorescent light source is then applied to the microarray to generate a fluorescent 
image identifying which elements of the microarray bonded with fragments of the DNA sample and which did not. The 
image is then analyzed to determine which specific DNA fragments were contained within the original sample and to 
determine therefrom whether particular diseases, mutations or other conditions are present within the DNA sample. 
[0005] For example, a particular element of the microarray may be provided with fragments of DNA representative 

30 of a particular type of cancer. If that element of the array fluoresces under fluorescent illumination, then the DNA of 
the sample contains the DNA sequence representative of that particular type of cancer. Hence, a conclusion can be 
drawn that the patient providing the sample either already has that particular type of cancer or is perhaps predisposed 
towards that cancer. As can be appreciated, by providing a wide variety of known DNA fragments on the microarray, 
the resulting fluorescent image can be analyzed to identify a wide range of conditions. 

35 [0006] Unfortunately, under conventional techniques, the step of analyzing the fluorescent pattem to determine the 
nature of any conditions characterized by the DNA is expensive, time consuming, and somewhat unreliable for all but 
a few particular conditions or diseases. FIG. 1 illustrates various conventional techniques for analyzing the flourescent 
pattem. Some prior art systems utilize two or more of the techniques. It should be noted though that many combinations 
of the techniques are not provided in the prior art. Briefly, the fluorescent pattem is quantized to yield a dot spectrogram 

40 at step 10. Any of four different techniques are employed to identify oligonucleotides represented by the dot spectrogram 
and to then identify mutations based upon the oligonucleotides. More specifically, the dot spectrogram may be analyzed 
using a trained neural networi< recognizer 12, a statistical decision theory recognizer 14, a fuzzy/expectation method 
(EM) clustering algorithm 16 or a rule-based inferencing/tmth table search 18. 

[0007] The results are interpreted to yield mutations of Interest at step 20. Then the mutations are again processed 
45 using either a trained neural network recognizer 22, a statistical decision theory recognizer 24, a fuzzy/EM clustering 
algorithm 26 or a rule-based inferencing/truth table search 28. The results are combined at step 30 to yield a diagnosis. 
[0008] Finally, disease confirmation is performed at step 32 by "reduction", in other words the disease is confirmed 
by probabilistic inferenctng. 

[0009] One major problem with many conventional techniques is that the techniques have poor repeatabil ity. Hence, 
50 if the same sample Is analyzed twice, different results are often obtained. Also, the results may vary from lab to lab. 
Also, skilled technicians are required to operate the DNA microarrays and to analyze the output resulting in high costs. 
One reason that repeatability is poor is that the signatures within the dot spectrogram that are representative of mu- 
tations of interest are typically very weak and are immersed in considerable noise. Conventional techniques are not 
particularty effective in extracting mutation signatures from dot spectrograms in low signal to noise circumstances. 
55 [0010] Accordingly, it would highly desirable to provide an improved method and apparatus for analyzing the output 
of the DNA microarray to more expediently, reliably, and inexpensively determine the presence of any conditions within 
the patient providing the DNA sample. It Is particularty desirable to provide a technique that can identify mutation 
signatures within dot spectrograms even In circumstance wherein the signal to noise ration is extremely tow. It is to 
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these ends that aspects of the invention are generally drawn, 

1001 1] One analysis technique for achieving the aforementioned advantages is described in U.S. Patent 6,1 36.541 
entitled "Method and Apparatus for Interpreting Hybridized Biochip Patterns Using Resonant Interactions Employing 
Quantum Expressor Functions". Briefly, the method of the co-pending application operates as follows. The method 

5 identifies mutations, ifany. present in a biological sample from a set of known mutations by analyzing a dot spectrogram 
representative of quantized hybridization activity of oligonucleotides in the biological sample to identify the mutations. 
A resonance pattem is generated which is representative of resonances between a stimulus pattern associated with 
the set of known mutations and the dot spectrogram. The resonance patlem is interpreted to yield a set of confirmed 
mutations by comparing resonances found therein with predetennined resonances expected for the selected set of 

10 mutations. In a particular example described in the co-pending application, the resonance pattem is generated by 
iteratively processing the dot spectrogram by perfonming a convergent revertaeration to yield a resonance pattem rep- 
resentative of resonances between a predetemnined set of selected Quantum Expressor Functions and the dot spec- 
trogram until a predeteimined degree of convergence is achieved between the resonances found in the resonance 
pattem and resonances expected for the set of mutations. The resonance pattem is then analyzed to a yield a set of 

15 confirmed mutations by mapping the confirmed mutations to known diseases associated with the pre-selected set of 
known mutations to identify diseases, if any, indicated by the DNA sample. A diagnostic confirmation is then made by 
taking the identified diseases and solving in reverse for the associated Quantum Expressor Functions and then com- 
paring those Quantum Expressor Functions with ones expected for the mutations associated with the identified disease 
to verify correspondence. If no correspondence is found, a new sub-set of known mutations are selected and the steps 

20 are repeated to detenmine whether any of the new set of mutations are present in the sample. 

[0012] By exploiting a resonant interaction, mutation signatures may be identified within a dot spectrogram even in 
circumstances involving low signal to noise ratios or. in some cases, negative signal to noise ratios. By pemnitting the 
mutation signatures to be identified in such circumstances, the reliability of dot spectrogram analysis is thereby greatly 
enhanced. With an increase in reliability, costs associated with perfomriing the analysis are decreased, in part, because 

25 there is less of a require ment for skilled technicians. Other advantages of the invention arise as well. 

[0013] Although the method of the co-pending application represents a significant advance over techniques of the 
prior art, room for further improvement remains. In particular, it would be desirable to enhance the method of the co- 
pending application to achieve a higher degree of repeatability. 

[0014] In particular, repeatability is affected by spatio-temporal degradation of hybridization in DNA microan-ays im- 
30 plementing both passive and active hybridization. In bioelectronic systems implementing passive hybridization, sources 
affecting repeatability of analysis results include: 

Stochastic variability in chemical kinetics 
Immobilized oligonucleotide damage during fabrication 
35 . Uneven kinetics during thermally facilitated fluidics reaction 

Post hybridization themnal degradation. Cun-ently biochips are "amplification limited". 

This Is in large part due to losses during high-tennperature hybridization downstream. During periods when the sample 
temperature changes from high to low or low to high, extraneous, undesirable reactions can occur that consume im- 
40 portant reagents and create unwanted and Interfering chemicals. Rapid transitions ensure that the sample spends a 
minimum of time at undesirable intennediate temperatures, so that the amplified DNA product has optimum fidelity 
and purity. So current methods rely on excessive amplification to compensate for these losses. 

Oligonucleotide entanglement 
45 - Environmental decoherence due to energy and radiatton 
Uneven fluidic catalysis 

Unstable fluorescence and chemiluminiscence maricer binding 
Spontaneous emissions 
Partial bindings 
50 . Anti-aliasing during readout and digitization 

[0015] Active hybridization is degraded by 

Capacitlve coupling between elements of the Immobilized matrix 
55 - Partial bindings due to current leakage and uneven conductance 
Ultrascale quantum squeeze effects 
Spontaneous emission 
Nonspecific oligonucleotide trapping 
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Chaotic relaxation across the array 

[0016] Hence, aspects of the present invention are directed, in part, to providing enhanced repeatability of biological 
sample analysis despite these factors. 

5 

Summary of the Invention 

[0017] In accordance with the invention, a method according to daim 1 is provided for identifying mutations, if any. 
present in a biological sample. The method operates to analyze an anrayed microdevlce output pattem generated using 

10 the sample to identify the mutations in the sample. In accordance with the method the output pattem is tessellated. A 
stimulus pattem associated with the set of known mutations is generated. A resonance pattem is then generated which 
is representative of resonances between the stimulus pattem and the tessellated output pattems. The resonance pat- 
tem is interpreted to yield a set of confrnmed mutations by comparing resonances found therein with predetenmined 
resonances expected for the selected set of mutations. 

15 [001 8] In an exemplary embodiment, the output pattem Is a dot spectrogram representative of quantized hybridization 
activity of oligonucleotides in a DNA sample. The stimulus pattem is generated based upon Quantum Expressor Func« 
tions. The dot spectrogram is tessellated to match morphological characteristics of the Quantum Expressor Functions 
and local parametrics are extracted. The tessellated dot spectrogram and the stimulus pattem are transfomned to a 
metrically transitive random field via phase shifting. The resonance pattem is generated by iteratlvely processing the 

20 tessellated dot spectrogram by performing a convergent reverberation to yield a resonance pattem representative of 
resonances between the Quantum Expressor Functions and the tessellated dot spectrogram until a predetermined 
degree of convergence is achieved between the resonances found in the resonance pattern and resonances expected 
for the set of mutations. The convergent reverberation includes the step of performing a convergent revert>erant dy- 
namics resonance analysis of the tessellated dot spectrogram using the resonance stimulus pattem to identify muta- 

25 tions represented by the tessellated dot spectrogram. The convergent reverberation also includes the step of performing 
a convergent reverberant dynamics resonance analysis of the mutations using the resonance stimulus pattem to iden- 
tify diagnostic conditions represented by the mutations. 

[001 9] Also in the exemplary embodiment, the convergent reverberant dynamics resonance analyses are performed 
by determining resonance dynamics relaxation values based upon the tessellated dot spectrogram and the resonance 

30 stimulus; filtering the dynamics relaxation values using ensemble boundary and complete spatial randomness (CSR) 
filters to yield a second set of values; applying bulk property estimators to the dynamics relaxation values to yield a 
third set of values; evaluating the second and third sets of values to determine a degree of resonance convergence; 
and then determining from the degree of resonance convergence whether a paralysis of dynamics has occun-ed and. 
if so. repeating the aforementioned steps. 

35 [0020] In the exemplary embodiments, by tessellating the dot spectrogram to match morphological characteristics 
of the Quantum Expressor Functions and by exploiting a resonant interaction employing a resonance convergence 
check which uses extracted tessellation parametric, mutation signatures may be identified within a dot spectrogram 
with a high degree of repeatability. By achieving a high degree of repeatability, the reliability of dot spectrogram analysis 
is thereby greatly enhanced. With an increase in reliability, costs associated with performing the analysis are decreased. 

40 in part, because there Is less of a requirement for skilled technicians. Other advantages of the invention arise as well. 
[0021] Among other applications, principles of the invention are applicable to the analysis of various arrayed biomo- 
lecular, Ionic, bioelectronic. biochemical, optoelectronic, radio frequency (RF) and electronic microdevices. Principles 
of the invention are particulariy applicable to mutation expression analysis at ultra-low concentrations using ultra-high 
density passive and/or active hybridization DNA-based microarrays. Techniques implemented in accordance with the 

45 invention are generally independent of the physical method employed to accumulate initial amplitude information from 
the bio-chip array, such as fluorescence labeling, charge clustering, phase shift Integration and tracer imaging. Also, 
principles of the invention are applicable to optical, optoelectronic, and electronic readout of hybridization amplitude 
pattems. Furthemnore. principles of the inventton are applicable to molecular expression analysis at all levels of ab- 
straction: namely DNA expression analysis. RNA expression analysis, protein interactions and protein - DNA interao- 

50 tions for medical diagnosis at the molecular level. 

[0022] Apparatus embodiment according to claim 28 is also provided. 

Brief Description of the Drawings 

55 [0023] The features, objects, and advantages of the present invention will become more apparent from the detailed 
description set forth below when taken in conjunction with the drawings in which like reference characters identify 
correspondingly throughout and wherein: 
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FIG. 1 is a flow chart illustrating exemplary techniques for analyzing dot spectrograms, some combinations of 
which are in the prior art. ^r^^,. 
FIG. 2 is a flow chart illustrating an exemplary method for analyzing the output of a hybridized DNA microarray 
biochip in accordance with the method of the «>i)ending application. 
5 FIG. 3 is a flow chart illustrating exemplary enhancements to the method of the co-pending application provided 

in accordance with the present invention. . 
FIG. 4 is a flow chart illustrating an exemplary method for analyzing the output of a hybndized DNA microarray in 
accordance with the method of the present invention. 

10 Detailed Description of Exemplary Embodiments 

100241 With reference to the remaining figures, exemplary embodiments of the method of the co^)ending applicafion 
and of the present invention will now be described. The exemplary methods will be described primarily wrth respect to 
the analysis of mutation signatures within output patterns of DNA biochip microarrays. but principles of the invention 
15 to the analysis of a wide variety of other patterns as well. 

Overview of Co-Pending Application 



20 



25 



30 



r0O251 Briefly, Bie exemplary method of the co-pending application exploits, among other features: (a) a novel rep- 
resentation, interpretation and mathematical model for the immobilized oligonucleotide hybridization patterns, repre- 
sented via a dot spectrogram: (b) a new "active" biomolecular target detection and discrimination method based on 
quantum resonance interferometry. and (c) a new spatial hashing function that yields accurate diagn^tic assessment. 
[00261 To this end the exemplary method of the co-pending application exploits a fundamentally different computa- 
tional paradigm for mutation expression detection in pre-enhanced dot spectrogram realizations. The method is an 
innovative modification to dynamically arrayed quantum stochastic resonance (QSR) for discrete system analysis. The 
arraying strategy is a function of the expression pathway of interest. The method depends on the molecular diagnostic 
spectmm being addressed. Banks of coupled quantum resonators are algorithmically designed to significantly enhan^ 
signal-to-noise (SNR) perfomiance and fuse multiple synthetic renomralized dot spectrogram realizations to better 
detect prespecified biomolecular expression patterns. 

100271 Moreover, the exemplary method of the co-pending application exploits an enhancement in previous exten- 
sions to classical stochastic resonance (SR) and anray enhanced SR (AESR) in signal Processing and censor da a 
analysis. Stochastic resonance is a phenomenon v*erein the response to a sensor, modeled in tenris of a bistable 
nonlinear dynamical system, is enhanced by applying a random noise element and a periodic sinusoidal forcing func- 
tion SR occurs wrhen the SNR passes through a maximum as the noise level is increased. 
35 [00281 Thus as important aspect of the exemplary method of the couending application involves the coupling of 
transfonned and preconditioned discrete microarray outputs to a mathematical model for a quantum-mechanical dy- 
namical system with specific properties. When driven in a particular manner, the coupled system exhibits a nonlinear 
response that corresponds to detection of phenomena of interest. The method exploits modulation of observables frorn 
a -base- (canonical continuous dynamical system), so that a selected set of spectral properties match a sinriilar selected 
40 spectral properties of a discrete spatial tessellation substmcture from an amplitude spectrogram derived frorri bioelec- 
tronic observables. The method further exploits the concept of convolving a discrete spatial system (denved from base 
mutants of interest) with a continuous asymmetric temporal system to derive a spatiotemporal input to further convolve 
with another discrete spatial projection (of an inherently partially stabilized spatiotemporal system). 
[00291 Hence key components of the exemplary biomolecular detection method of the co-pending are: (i) selection 
45 of a basis system: (ii) generation of designer Quantum Expressor Function (QEF) for coupling with the substrate to be 
analyzed- (iii) generation of a Hamiltonian to describe relaxation dynamics of the coupled system; (iv) tnodulation of 
reswiance parameters to enforce earty resonance: (v) and exploitation of resonance suppressors to verify detection. 
(00301 Referring to FIG. 2. initially at step 1 00. a set of mutations of interest are selected. The mutations, for example, 
may be mutations relevant to cancer. AIDS, or other diseases or conditions. At step 101. preconditioner transfomis 
50 are oenerated based upon the selected set of mutations. The preconditioner transfomis are provided to convert mu- 
tation nucleotide sequences into expected amplHude patterns in the prespecified microarray representation, given a 
particular biochip layout. At step 102. Quantum Expressor Functions are generated based upon the Hamiltonian of a 
pre-selected basis system. The Quantum Expressor Functions are designed to couple the Hamiltonian for the selected 
basis system to a predetennined DNA microarray configuration to pemiit a resonance interaction involving the output 
55 of the DNA microan^y. Resonance stimulus is generated, at step 106, using the Quantum Expressor functions. 

[00311 What has been summarized thus far are preliminary steps performed off-line for setting up the Quantum 
Expressor Functions and the corresponding resonance stimulus. These steps need be perfomied only once for a given 
set of mutations and for a given DNA microanay configuration. Thereafter, any number of output patterns from the 
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DNA microarray may be processed using the Quantum Expressor Functions to identify whether any of the mutations 
of the pre-selected set of mutations are found therein. Preferably. Quantum Expressor Functions pre^enerated 
for a la%e set of mutations and for a large set of DNA micnjarray patterns such that, for each new DNA ni'«-oarray 
output ^ttem from each new patient sample, the presence of any of the mutations can be quicWy identified using the 
predetermined set of Quantum Expressor Functions. In general, the aforementioned steps need 
update the Quantum Expressor Functions to accommodate new and different DNA microarray patterns or to if new 
mutations of interest need to be considered. ... _!-^«,rov, 

[00321 At step 1 06. an output pattem (ref en^d to herein as a Dot Spectrogram) is generated using a DNA microarray 
or which Quantum Expressor Functions have already been generated. At step 108. the dot spectrogram is precond^- 
tioned to yield a dot spectrogram tesselation (DST) to permit exploitation of a resonance between the ^t s^ctrogram 
and the Quantum Expressor Functions. The actual resonant interaction, which involves convergent raverberatwns. is 
perfomied at step 110 until a pre^etem^ined degree of convergence is achieved. 

resulting resonance pattem is processed at step 11 2 to identify any mutations represented theraby. wiH be desonbed 
below step 112 is rendered trivial by virtue of the aforementioned resonant interaction which is based upon Quanhim 
Expt^ Function already correlated with the pre-selected mutations. Hence, no complicated ana^rsis is required to 
interpret the resonance pattem to identify the mutations. Next, at step 114. the mutations are mapped to corresponding 
diseases and conditions to thereby identify any diseases or conditions that the patient providing the sample being 
analyzed is afflicted with. Again, this is a fairly trivial step. Finally, at step 116. diagnost-c confinriaUon pre«o^ed to 
verifythat the diseases or conditions are present in the sample. This is achieved by starting with the found d«eases 
or conditions and then perfomiing the steps of the method in reverse. „^ H»f=iiorf 

[0033] Each of the aforementioned steps are described in detail in the co-pending application and the detailed de- 
scription thereof is not repeated herein. 

Overview of the Exemplary Method of the Presen t Invention 

[0034] The present invention is directed, in part, to improving the repeatability of the method of the co-p^ing ap- 
plication by tessellating the dot spectrogram so as to match morphological characteristics of the Quanhim Expr^or 
Functions and by using extracted local parametrics as part of a resonance convergence check dunng the ^sonan^ 
interaction. These additional steps have the advantage of establishing uncertainty bounds which permit method re- 
peatability to l)e enhanced and quantified. 

[0035] FIG. 3 niustrates enhancements to the technique of FIG. 2 provided by the present invention along with steps 
of the technique of FIG. 2. The repeated steps of FIG. 2 appearing in FIG. 3 may be the same as 't^o^e of FIG. 2 and 
will not be redescribed. Like reference numerals, incremented by one hundred, are employed to represent the repeated 
steps. 

[0036] Briefly the enhancements of the present invention entail five major steps: 

. Preconditioning the hybridized array output pattem (i.e. the dot spectrogram) by a fuzzy tessellation and coupling 
the preconditioned output pattem with a canonical system with aftereffect and memory properties (step 218); 

- Estimation of amplitude wanderings or dispersion (step 220); 

- Implementation of resonance interaction by integrating partial or subthreshold resonances using phased an^y 
enhancement operator resonance dynamics (step 224) with additional resonances synthetically induced to ac- 
commodate the possibility for the presence of single-point and two-point mutations around the mutation-centered 
pixels; and 

Resonance recombination (step 226). 

[0037] A combination of one or more these enhancements are superimposed on the techniques described in the co- 
pending patent application to address specific sources of hybridization degradation, device imperfections and protocol 
variability in the analysis process to thereby enhance repeatability. 

[0038] Referring now to FIG. 4. initially at step 300. a set of mutations of interest are selected and at step 301 
preconditioner transfomis are generated based upon the selected set of mutations. At step ^O^. Q^-turn E^^^^^^^^ 
Functions are generated based upon the Hamiltonian of a pre^Iected basis system. Phase shifted resonance stimulus 
is oenerated at step 304. using the Quantum Expressor functions. Grouping stimulus is also generated, at step 305. 
100391 Steps 300 - 305 are preferably perfomied off-line to set up the Quantum Expressor Functions and the corre- 
Und ng resonance stimulus and grouping stimulus and need not be repeated other than to "P«lat« Q^f"'""" 
Btpres^r Functions to accommodate new and different DNA microarray patterns or if new mutations of interest need 

!oMOr"At'S*306, a dot spectrogram is generated using a DNA microanay for which Quantum Expressor Funrtions 
have already been generated. At step 307. the dot spectrogram is tessellated to match morphological characteristics 
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of the Quantum Expressor Functions yielding a dot spectrogram tesselation (DST). At step 308, local parametncs of 
the tessellated image are extracted. Then, at step 31 0. an amount of amplitude wandering is detemiined and compared 
with pr&<letermined allowable generator function limits. If. at step 310, the amplitude wandenng is not within the al- 
lowable generatorfunction limits, then execution proceeds to step 312 where the tessellated dot spectrogram to match 

s spectral characteristics of the Quantum Expressor Functions. Steps 308 and 310 are repeated un il the amplitude 
wandering is found to be within the pre-determined limits at which point execution proceeds to block 314 wherein a 
resonance interaction is perfoimed between the tessellated, renomnalized dot spectrogram and the phase shifted res- 
onance stimulus generated at step 304 and the group stimulus generated at step 305 to identify any mutations repre- 
sented by the dot spectrogram. . ^ ^ . „ ■ ^ 

10 [00411 The actual resonant interaction, which involves convergent fBverberattons. includes the following sub-steps 
also shown in FIG. 4. At step 316, a resonance dynamics iteration is initiated v*ich includes the use of ensemble 
boundary and CSR operators (step 318) and the use of bulk property estimators (step 320). The ensemble boundary 
filters CSR filters and bulk property estimators are applied to the tessellated, re-nomialized dot spectrogram in com- 
bination with the resonance and group stimulus. The resulting filtered dot spectrogram is then evaluated to delemiine 

15 a deqree of resonance convergence to one or more of the set of predetemiined mutations, at step 322. The degree of 
convergence is evaluated, at step 324, against a Lindbald condition and. if the Lindbald condition Is not met. the system 
is deemed to be subject to paralysis of dynamics and execution proceeds to step 326 wherein possible hixel death is 
compensated for by increasing a time scale for the iteration initiated at step 316 and then repeating the iteration. 
10042] Here it should be noted that mutation death and paralysis of dynamics are different concepts. The mutation 

20 death check is a condHtonal check . If this check shows that a resonance is not possible for a specific mutation reso- 
nance centered (MRC)-hixel then the iteration is terminated and block 314 is exited. But failure of resonance dynanriics 
is not sufficient to conclude that a specific mutation is absent. Indeed, if the "hixel death" check fails, that implies that 
resonance could be still obtained in a downstream iteration. 

[00431 If at step 324. the iteration has converged and no paralysis has occurred, then one or more mutatioris have 
25 been reliably Identified. Execution proceeds to block 325 wherein another resonant interaction is perfomied to identity • 
particular diseases represented by the mutations. Briefly, at step 326. a resonance dynamics iteratfon is initiated which 
includes the use of ensemble boundary and CSR operators (step 328) and the use of bulk property estimators (step 
330). The resulting filtered dot spectrogram is then evaluated to determine a degree of resonance convergence, at 
step 332. The degree of convergence is evaluated, at step 334, against a Lindbald condition and. if the LindbaW 
condition is not met. the system is deemed to be subject to paralysis of dynamics and execution proceeds to step 336 
wherein a time scale for the Iteration initiated at step 326 is Increased before the Iteration Is repeated. If the Lindbald 
conditton is met. then diseases corresponding to the mutations found using step 314 have been reliably identified. 
[00441 Processing in accordance with step 335 depends on the blochip. The flowchart of FIG. 4 illustrates a general 
fonn of the method requiring processing during step 335. In other implementations, this step is trivial or can be elimi- 
nated entirely. In this regard, the overall method is implemented at two levels of abstraction, depending on how well 
the disease genomics is understood. Detection of a specific mutation is necessary and sufficient to conclude expressiori 
of a specific gene. But the expressed gene may not be the one to conclusively identify the disease. Then another level 
of abstraction is invoked wherein the method is applied inferentially by expanding the gene expression circuit or gene 
expression tree to determine if there is evidence that all expressed genes eventually lead to one that conclusively 
<o Identifies the disease. So step 335 operates on the results of step 314, such that all identified mutations are used ^ 
an input to detennine if the complete expression pathway for leading up to the point that a disease can be concluded. 
[00451 If the biochip is so designed that the mutations con-esponding to all intemiediate expressed products, from 
any disease starting point can be captured by resonance output of step 314. then sub-steps within step 335 and sub- 
sequent steps 340 and 342 can be circumvented. If not. clustering step 340 and geometnc hashing step 342 are 
45 provkled to identify that an expression pathway is present that trivializes the disease conclusion =(step 344). 

[00461 Once the diseases are identified, clustering properties are evaluated at step 340 to selectively eliminate oli- 
aonudeotWes representing possible diagnoses based on morphological filtering of subthreshold resonances and any 
subsequent recentering (i.e. the inverse of dispersion). Steps 314-335 produce a cluster of sub-threshold resonances. 
Step 340 is a reveriftaation such that all induced resonances are present in the target sample and not a manifestation 
so of multiple rescaiings and synthetic SNR enhancements. ....... j ^ 

[0047] Then a hashing projector is applied at step 342 to order the mutations. A diagnostic decision is then rendered 
at step 344 by examining the order of the mutatkws and comparing the mutattons with a table identifying con-esponding 
diseases 

r00481 Thus the output of block 314 represents all hixels that identify complementary oligonucleotide bindings in the 
55 biological sample being analyzed and this represent "mutations". The output of block 335 comprises a set of expressed 
genes that are associated with a particular pathogenic pathway and thus represent a preliminary "diagnosis . Further 
analysis of the pathogenic pathway provides a set of possible diseases, if any. This decomposition is motivated by 
scaling the computation to answer three questtons: 



30 



35 
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- Does the current condition imply predisposition to a specific disease? 

: Ta! re r ^^SSitr; ^.e o^^..^ fr^ .ne targe. san,p,e. given a specfc genomic 
encoding implemented by the biochip? 

100491 In anv case if the diagnostic decision rendered at step 344 is afRm^ative. then the diagnosis is output. If the 
LonLs^nSre thrn at sLp 346. a determination i^ 

whrt:rn«i:?:rofm'^^^^^^^^^ 

made by examining a table listing all possible mutations. If there are aHematve mutat«ns. then '^^'^ 
from step 300 If not. then a signal is simply output indicating that no mutations were found in the sarnple. 
S NOW dSs of the steps of the new method will be provided. Details regard steps already descnbed .n the 
co-pending application will not be repeated herein. 



Mutations Sets 



10051] The mutation set of interest generated at step 300 is selected by identifying ol-gonuclaot.des represen^fcve 
S the {Z} mutations of interest. Each oligonucleotide is represented by which .s g«/en by " 
W. C T Gl base associated with each array cell |a.b] where 10<kS25. The entire set of unique oligonucleotides denoting 
I^SJS of int'STS). is given by where 8 = {A.C.T.G, length |A|= 'ength M. and 0<|^ - ,;^| . k. and 

the desioned in vi/(l m) oligonucleotide sequence is a perfect complement to only A(l,m) for all l.m. 

00521 AfpLrtof^ep 300, an oligonucleotide table is generated which contains the oligonudeoWe sequences as- 

ocSldCth^a^ mutation Of InterLt identified by r«w and column 1^^^^^ 

for subsequent use at step 312 to map locations within the dot spectrogram ^XTwso^oaJS sS? 

oligonucleotides such that mutations present in a sample being analyzed are easily ^ 
a mutation table is generated which contains the diseases associated with each mutation of interest. The mutahon 
faWe fs provSed for subsequent use at step 314 to map mutations identified at step 312 to specific diseases or other 
medical conditions such that the diseases can be easily identified. 

Basis System for Quantum Expressor Functions 

[0053] The selection of the basis system and the generation of the QEF's based thereon depends Pf • f 
Siaracteristicsof the DNAmicroarray. Inthe exemplary embodiment, the DNAmicoan-ay is an N by M DNA chip array 
wherein an element of the array is refened to herein as an ^toxeT; o(iJ). 
[0054] The pre-hybridization microarray (PEBC) is expressed as: 



■ I 



where N and M refer to the linear (row and column) dimensions of the 2-D microarray. 
[0055] The numeric value associated with each oxel is given by: 

<5(/. j) = a/r^*"' +tt*.i •4*"^+...+ai .4^+00-4° 

where [al=IA|C|T|Gl take the values [0|1 12131 respectively 
100561 An element of the dot spectrogram is refened to herein ^ a hixel: h(ij). 

!S rsSnSnbasissystemisLectedforusewiththistypeofarray.Otherbasissystemmaybeappre^^^^ 
for either the same or other microan-ay configurations.. 

Quantum Expressor Functions 

10058] The QEF is generated at step 302 based upon the spm Boson basis system by first "f,-:"!"^^^^^ 

orTesystem. calculating hamionicamplitudeslPJfor the Hamiltonian.generatng an <^^^^ 

eSnment e^tramment states of the OF of the greund truth and finally modulating the OF of ground tr«th to yield the 

fo059] The QEF-s generated at step 302 are converted to a phase-space representation. Also, if the output of the 
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hybridization chip is not in phase space then it is converted as well. The conversion ,s perfonr-ed us ng Phase emb^ 
ding operator, r. described in the copending application. Results associated with comb.natonal H°P' Algebra are used 
to Mrtain amplitude dispersion due to loss of hybridization. A special case of quantum ^^''V.^f f *':'!!'1^„^ 
Segal (GNS) construction is used to disperse group stimulus. Note that coproduct construct of the Hopf algebra plays 
tt^role of "sharing our possible explanations of a fact. The GNS dispersion of QEF is implemented using an approx- 
imation: 

<t>Q£p{ainp.vector) = 0,30]' where 0, = e"*" 

where H denotes the hamiltonian for the coupling spin boson system. 
Generation and Tesselation of the Dot Spectrogram 

[00601 As noted a dot spectrogram is generated at step 306 for a sample from an N by M DNA chip array wherein 
an element of the array is an "oxe/"; o(U)- A 6-c manufacturing process accuracy in microarray design is assumed. 
Each array cell amplitude is giv"i;ib7«I>(i.j) for i: 1 to N. and j: 1 to M. Let denote the a pnon known oligonucleotide 
given by [ooOvOj. where a = {A.C.T.G} base associated with each an^y cell [a.b] where 1 0<k<25. The complimentary 
strand, derived fix)m unknown sample is denoted by 7). « » 

100611 The post-hybridization microarray is treated mathematically using the machinery of equations with afterelfed 
Each hixel given by *(i.j) is represented as a cluster of dynamical systems of potentially ICB] correctly bound. [UB] 
unbound, (PB] partially bound and [IB] incoffectly bound. Thus [CBl^j, + WJBlt^n) * [PB]*oj) * tlB]*(ij) = Vm) 

[0062] ^°The dot spectrogram *(i.j) is then tessellated to determine idealized ensemble boundaries for forcing down- 
stream resonant action. . . . . , . . ,„ 

[00631 Typically, in signal processing applications, high pass or band pass spatial filtenng is implemented to enhance 
SNR in DS matrix. Altemate methods apply a combination of Laplacian or other edge detection filters appfy to enhance 
signal from arrays cells with a higher hybridization concentration from those of the adjacent cells. These SNR enhance- 
ment methods however work only with positive or zero-SNR. Since SNR in general is negative in our case (ulfra-low 
target DNA concentrations), these methods in effect amplify noise or further blur the hixel boundanes. 
[00641 Tesselation is performed by performing gradient refocusing and reseating as described in the oo-pending 
applicatton. In the altemative. a Dirichlet tessellation operator or a Delaunay triangulation operator are applied to tes- 
sellate the dot spectrogram. 

Extraction of Local Parametrics and Calculation of Amplitu de Wandering 

[0065] The tessellated image is treated as a metrically transitive random field. All properties associated with a singular 

(detemiinistic), homogeneous (i.e., stationary) field are subsumed. 

(0066] The parametric of most interest is the integrated density of states, given by 



where 



whence 
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where n is the number of eigenvalues to the system (random field approximation). 
[0067] This is computed for each tessellation region in the dot spectrogram. 

[0068] Amplitude wandering is determined using Palm generators as described in the co-pending patent #1. The 
Palm generators exploits the notion of generator functions to capture stochastic variability in hybridization binding 
efficacy. The exemplary method described herein draws upon results in stochastic integral geometry and geometnc 
probability theory. . 
[0069] -Amplitude wandering estimate" that bounds the hlxel amplitude dispersion due to total hybndization losses, 
is computed using Palm generators over the globally re-scaled dot spectrogram to capture amplitude wanderings and 
transitions at element, neighboring pair and local ensemble levels. Step 310 provides a measure for each mutation- 
recognizer centered (MRC-) hixel that is invariant to local degradation. The measure Is expressed via the fbmn 

where Z denotes the set of mutations of interest. In other words, we determine the function f(z) under the condition 
that m(z) should be invariant with respect to all dispersions ^. Also, up to a constant factor, this measure is the only 
one which is invariant under a group of motions in a plane. In principle, we derive detemninistic analytical transfonnations 
on each MRC-hixel., that map en-or-elliptic dispersion bound defined on^R^ (the two dimension Euclidean space - i.e.. 
oxel layout) onto measures defined onSR • The dispersion bound is given by the fornn 

Log4(6(,-^|^). 

[0070] Recall that Palm distribution, n of a translation (T^) invariant, finite intensity point process in«R" is defined to 
the conditional distribution of the process. It is expressed in terms of a Lebesgue factorization: 



[0071] Where n and A completely and uniquely detemnine the source distribution P of the translation invariant point 
process. The term Ep N* denotes the first moment measure of the point process and is a probability measure. In 
patent #1 we described how to compute n and A which can uniquely encode the dispersion and amplitude wandering 
associated with the MRC-hixel. 

[0072] In this invention we relax the strong assumption that Palm generators, n and A. capture ail sources of sto- 
chasticity in dot spectrogram output. Since hybridization losses are affected in unknown and unpredictable ways, we 
need to modify the generators as probabilistic functions themselves. In other words the generators are converted to 
manifolds as opposed to a point function. 

[0073] {Pmiijy <^m^^m*^m) specifies a continuous probability density function for amplitude wandering in the m-th 
MRC-hixel of interest where the tenms denote: 

oligonucleotide density per oxel p^,j), PCR amplification protocol (a^). fluorescence binding efficiency (ti J and 
imaging perfonnance ( ra^]). Previously we required a preset binding dispersion limit to be apriori provided to com- 
pute A, given by the second moment to the function at SNR = 0. 

[0074] Nondispersive n is computed using f! = 0 * P where 

[0075] And and represent nomialized hybridization dispersion limits (typically preset to 0.1 and 0.7 respectively 
to assume losses between 10% - 70% hybridization. 

[0076] Precond itioned dot spectrogram is represented by o (/,y)- where function 1 /(1 +exp( p (...))) was used to express 
the underiying known and stationary point process. 

[0077] The latter assumption is relaxed in this method and detenmination of whether the amplitude wandering is 
within allowable generator function limits is achieved by: 
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where 

Hn denotes the nondispersive generator. . ^ ^ 

^1 ^2. . provide the Laplace characteristic functional of the Poisson random field associated with each source 
of hybridization degradation. The contributions are estimated using 

7|, = Cd.(det5/)"^Xo ^ 

where c is gain constant and a denotes a nonrandom matrix 

§ = {afA<iJ<d) 

given by^ ^^^^^^^^^^^^j metrically transitive fields representing the unique solution of the following variational 

problems: 



(i) 



^' V/c 

[00781 The differential operator for the metrically transitive field for convoling the uncertainty parameters is denoted 
by Ao. Population and solution of the above equation requires esfimates for the fonA^ard sensitivity matnx of vanables 

impacting hybridization degradation. * x. -*u *.«^/.»roi 

[00791 Renonnalizatlon at step 312, if necessary, is perfonned on the tessellated .mage to further match spectral 
properties of the stimulus pattern The re-nomialization of the dot spectrogram is achieved by resealing the dot spec- 
trogram in the interval [-n. +k]. The entire calculation proceeds in the phase space which is why we transfonmed the 
system to the metrically transitive random field. 

Resonant Interaction to Identify Mutations 

[00801 As noted, at step 314. the resonant interaction between the QEF and the tessellated, re-nomialized dot spec- 
trogram is perfonned unfil a pre-selected degree of convergence is achieved. Resonance dynamics relaxation values 
are calculated at step 31 6 as follows. 

[0081] A closed-form convolutionless evolution equation is given by: 



where 
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DST Dsr, 



where both depend upon the nonnalized DST, (i.e.. initial state) at time t - post-hybridization but pre-conditioned 
state. And i'(f.t) is the unitary evolutere-«-'>« 

for the hamittonian specifying the after affect basis system dynamics. 
[0082] Also. 



10 



15 



20 



Dsr, S iDsr J 



[0083] In practice a small e.typically 1 0-6 is used. 

[0084] So if theoretical convergence time is Tq (outer convergence cycle time) and choosing t > x + Xq, then: 



^ (/) =6>(/,r+ro) 



25 

and 



30 



6>(r.r')=-t[^j[9(t>j 



and X is a time-ordering operator. 
35 [0085] The dynamics relaxation values are then filtered at step 31 8 using ensemble boundary and CSR filters (higher 
order Poisson kernel) as follows: 



40 



DST, 



where 



45 Pfie\)=—^ 2 

1-2rCos+r'^ 

where r^O. 

[0086] The bulk property estimators of step 320 are applied to the dynamics relaxation values as follows: 



50 



where t is the discretization step. 
[0087] The above expression provides an estimate of when a geomteric motion embodied by the convolutionless 
equation, is no longer a plausible resonance candidate. 

[0088] This is the closed form for an expression at which the coupling between DST and the microarray is broken 
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and a coupling with a nonlinear infonnation filter (NIF) is established. In essence, the system forgets any initial corre- 
lation and tends to a Lindbald condition. 

[0089] The resonance convergence Is determined at step 322 as follows: 



Log 



|U0-^1)-UjX[/^VGll 



> 1 



[0090] The system oscillates if no convergence is reached. If increasing the timescale x-times (-5) does not meet 

10 the condition, then the mutation is deemed to be absent. 

[0091 1 It should be noted that, unlike the technique of the co-pending application, in the present invention the absence 
of resonance over a maximum interation count does not imply absence of resonance. The reason is that both the dot 
spectrogram and the QEF are dispersed, i.e.. the SNR is reduced over an individual hixel, but is in fact increased over 
an ensemble. So the convergence decisions are made by cascading the inner loop reverberations as opposed to a 

15 single reverberation. So two timecycles are used for the convergence analysis: 

(a) time cycle over which hyperfine resonances are tracked, detected and used as a decision mechanisms to 
continue or stop the interation; 

(b) time cyde over which the absense of mutation is actually concluded. This is done by implementing a local 
20 maxima over output of previous step and then reintegrating.The method essentially accumulates partial resonanc- 
es and then applies the same resonance equation to the reseated and renormalized partial stage. 

[0092] This process can be analytically be represented as: 



30 



35 



40 



45 



50 



55 



> 0 



otherwise 



where c^ C2 and C3 are thresholding constants that are used to detect subthreshold resonances. Also, c^ » Cj » 
C3 » l/[amplitude resolution]. 

[0093] Also Ti and refer to the inner and outer integration timescales. In an implementation they refer to the iteration 
confer at which the integration loop is tenninated. exceeded or exited. Typically termination counter is set to one thou- 
sand steps with timescale of the order of ten nanoseconds for inner step and microsecond for outer step. So effective 
device convergence time is within onen hundred milliseconds for the entire computation. 

[0094] In this regard, if the Lindbald condition is not achieved and verified the dynamics is considered paralyzed. 
[0095] The paralysis of dynamics actually implies that coupling between 

(a) the Quantum Expressor Function, 

(b) the tessellated and nonnalized dot spectrogram , and 

(c) the convolutionless carrier 

is too weak to exhibit a nonlinear resonance. The physical interpretation is that the coupled system exhibits "fmstrated 
dynamics" which enhances and impedes resonance reaction at the same time. So the actual output takes the form of 
white noise over several hixels which oscillates. 

[0096] The detection of oscillation occurs when the spectral radius for the convergence criteria oscillates between 
limits [81 . £2] and does not tend towards O.This may be verified by tracking the spectral radius zero crossing with respect 
to the lower bound Zy If the zero crossing frequency exceeds a present number (e.g., 10) in this implementation, the 
dynamics is deemed paralyzed. 

If a paralysis of dynamics has occurred, a "mutation death" is evaluated as follows. The check for MRC hixel death 
relates to the verification of a suprathreshold resonance, where the resonance is defined as the integrand of partial 
resonances over the entire DST structure, i.e. 
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osr I 
> K 

V X1.T2 ^ predefined upper limit. Typically set to 100 for outer iteration and 1000 for inner iterations. 
[0097] The time scale for realization of the Lindt>ald condition is changed and the system reiterated. 
[0098] Hence the final output of step 314 is all hixels that identify complementary oligonucleotide bindings in the 
biological sample which are represented computationally by the set {hjj(i,j)} where {h|^(i,j)} is the corresponding oligo- 
nucleotide sequence looai . o^l for the kth surviving hixel. 

Resonant Interaction to Identify Diseases 

[0099] The mutations identified using block 314 are processed using similar steps within block 325 to identify diseases 
represented by the mutations. Hence the final output of step 335 is a set of expressed genes that are associated with 
a particular pathogenic pathway which is represented computationally by the set {VkC'j)) where {Vk(«.i)}- [aoai-aK] 
for the l-th element of the pathway capturing the k-th disease. 
[0100] For single disease analysis steps 324 - 334, i.e.. block 335 can be omitted. 

Clustering Property Check 

[0101] Diseases identified using block 325 are processed at step 340 to identify clustering properties as follows. The 
clustering operation is essentially a pruning operation based on morphological filtering of subthreshold resonances 
and subsequent recentering(i.e. the inverse of dispersion). 

[0102] The clustering computation is based on transversal ordering (is based on transversa! numbers) of the oligo- 
nucleotide sequencing undertying the resonance-centers for all subthreshold resonances. The concept draws from a 
result in hypergraph theory. Recall that transversal of a hypergraph H = {X: E^, Ej, . .. E^) is defined to be a set T <= 
Xsuch that T n E, 5t<> for I = 1. 2 m, where E^. E2, E^, define subgraphs. 

[0103] In this method, each oligonucleotide, associated with a mutation that survives "hixel-death" during resonant 
reverberation iterations, is represented by [ooav-otk]. where a = {A,C.T.G} base associated is treated as a sub- 
graph of the total set of unknown mutations that are actually present in the target sample. If the surviving hixel is an 
ensemble than each ensemble is treated as a subgraph with multiple nodes and several edges. If only an individual 
hixel survives than it Is treated as a single node subgraph. Transversal number of a hypergraph, H. is defined as the 
minimum number of vertices in a transversal. It is given by : 

f(iy) = mtn|r| 

[0104] The sub-steps involved in clustering are: 

- determine Ming ={Ai^2- »^^}- where A^. Aj. .., A,c denote the surviving resonance clusters. 
Next determine the following families: 



3,= 3.wU} -> rr{3,}= Min(rr3, v7V{^}) 



3,= 3,i^U} rr{3,}= Min(rr3, vrrU}) 



0 

otherwise 
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3^ =s — 

5 If Min A has k members, then the algorithm constructs Tr A = Tr \ in k steps. 

Hashing Projector 

[0105] A hashing projector is then applied at step 342 to the output oflhe clustering check. The hashing projector 
10 produces an enumeration of the leading k oligonucleotides with the highest transveral numl>ers. So a set of mutations 
or the corresponding expressed genes are created that have the highest sorted transversal numbers. Typically, all 
members that are seperated by a distance of, at most two, are chosen. 

Diagnostic Decision 

15 

[01 06] A diagnostic decision is rendered at step 344 based upon the output of the hashing projector. The diagnostic 
decision is achieved using a simple table lookup that is indexed by the results of hashing projection computation using 
the aforementioned tables. 

20 Alternatives 

[0107] Alternative possible mutations are evaluated at step 346. 

If alternatives are available, the alternative set of mutations of interest are loaded and the process is repeated beginning 
at step 300. Hence, if the original set of mutations from which the original set of QIF's were generated during the off- 
25 line process of steps 301 and 302. did not include the alternative mutations, then the off-line process is repeated with 
the new set of mutations to generate new QIF's. 

[01 08] In the event that method yields (and it often does) multiple disease detection hypotheses, all possible hypoth- 
eses are provided as plausible candidates. 

[0109] The technique described with respect to FIG. 4 is particulariy powerful in that it provides an enumerative 
30 solution which generally covers ait possible diagnostic candidates as opposed to only one or two, given the best ge- 
nomic understanding or mapping between expressed genes and diseases. 

Alternative Embodiments 

35 [01 1 0] Details regarding an implementation directed to measuring viral loads may be found in U.S. Patent 6,245,51 1 , 
entitled "Exponentially Convergent Therapy Effectiveness Monitoring Using Viral Load Measurements". 
[01 1 1 1 The exemplary embodiments have been primarily described with reference to flow charts illustrating pertinent 
features ofthe embodiments. Each method step also represents a hardware or software component for performing the 
corresponding step. These components are also referred to herein as a -means for" performing the step. It should be 

40 appreciated that not all components of a complete implementation of a practical system are necessarily illustrated or 
described in detail. Rather, only those components necessary for a thorough understanding of the invention have been 
illustrated and described in detail. Actual implementations may contain more components or, depending upon the 
implementation, may contain fewer components. 

45 

Claims 

1. A method for analyzing an output pattem of an arrayed microdevice to identify mutations, if any. present in a 
biological sample applied to the arrayed microdevice. said method comprising the steps of: 

50 

tesseilating the output pattem; 

generating a stimulus pattern associated with the set of known mutations; 

generating a resonance pattem representative of resonances between the stimulus pattem and the tessellated 
output pattem; 

55 interpreting the resonance pattem to yield a set of confimned mutations by comparing resonances found therein 

with predetemiined resonances expected for the selected set of mutations. 

2. The method of daim 1 wherein the output pattem is derived from an an^yed microdevice having a plurality of cells 
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each having a set of substantially identical immobilized oligonucleotides with the oligonucleotides of each cell 
being substantially unique from the oligonucleotides of every other cell. 

The method of daim 1 wherein said step of tessellating the output pattern comprises the step of amplifying highly 
local morphological variations in the output pattern. 

4 The method of daim 3 wherein the stimulus pattern is generated based upon Quantum Expressor Functions and 
wherein said step of tessellating the output pattern is perfbnned to tessellate the output pattem to match morpho- 
logical characteristics of the Quantum Expressor Functions. 

5. The method of daim 1 further induding the steps of: 

extracting local parametrics from the tessellated output pattem; 

detemiining whether a degree of amplitude wandering representative of the local parametrics is within a pre- 
determined allowable generator function limit; and 

if not, renonnalizing the tessellated output pattem to further match spectral properties of the resonance pattem. 

6. The method of claim 5 wherein said step of extracting local parametrics from the tessellated output pattem com- 
prises the steps of , ..^ XX I *• 

extracting parameters representative of an integrated density of states within each of a plurality of tesselation 
regions within the tessellated output pattem. 

7. The method of daim 5 wherein a degree of amplitude wandering is detennined by applying a Palm distribution to 
develop generators for estimating stochastic wandering. 

8. The method of daim I further induding the step of transfomiing the tessellated output pattem and the stimulus 
pattern to a metrically transitive random field. 

9. The method of daim 8 further including the step of renonmalizing the tessellated output pattern following transfor- 
mation to the metrically transitive random field. 

10. The method of claim 9 wherein said step of renormalizing is perfomfied to rescale the tessellated output pattem 
to the interval [-n, +n]. 

1 1 . The method of daim 1 wherein the step of generating a stimulus pattem assodated with the set of known mutations 
comprises the steps of: 

selecting a sub-set of the mutations for analysis; and 

selecting a sub-set of nonlinear Quantum Expressor Functions from a set of predetermined nonlinear Quantum 
Expressor Functions based upon the selected sub-set of mutations. 

12. The method of claim 11 further including the steps of: 

transfonning the nonlinear Quantum Expressor Functions and the tessellated output pattem into phase space. 

13. The method of daim 1 wherein the step of generating a resonance pattem indudes the step of 

iteratively processing the output pattem by perfomntng a convergent reverberation to yield a resonance pat- 
tem representative of resonances between a predetemiined set of selected Quantum Expressor Functions and 
the tessellated output pattem until a predetermined degree of convergence is achieved between the resonances 
found in the resonance pattem and resonances expected for the set of mutations. 

14. The method of daim 13 wherein the step of iteratively processing the output pattem by performing a convergent 
revertjeration indudes the step of perfonming a convergent reverberant dynamics resonance analysis of the tes- 
sellated output pattem using the resonance stimulus pattem to identify mutations represented by the tessellated 
output pattem. 

15. The method of daim 14 wherein the step of perfonming a convergent revert)erant dynamics resonance analysis 
indudes the steps of: . 
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a) determining a first set of resonance dynamics relaxation values based upon the preconditioned output 

^rri^lhrd;^"^.:.:^^:^^^^^ usmg ense^^e t^undary and CSR f.I.ers to yield a second set of 

c) aS'ying bulk property estimators to the dynamics relaxation values to yield a third set of values: 

d evaluates the second and third sets of values to determine a degree of resonance "'"^^^f 

e) determining from the degree of resonance convergence whether a paralyse of dynamics has occurred and. 

if so, repeating steps a) - e). 

16. The method of claim 1 5 wherein said step of determining resonance dynamics relaxation values based upon the 
preconditioned output pattern and the resonance stimulus comprises the steps of 

applying a convolutionless evolution operator to the tessellated output pattern. 

17. The method of daim 15 wherein said step of applying bulk property estimators to the dynamics relaxation values 
to yield a third set of values comprises the steps of „,.i,„,„„„«i 

applying a coupling operator to which couples the dynamics relaxation values representat«,e of the tessel- 
lated output pattern to a nonlinear information filter. 

1 8. The method of daim 1 5 wherein said step of evaluating the second and third sets of values to determine a degree 
of resonance convergence comprises the steps of ^. . ^ a«or =. nrpripfpr- 

determining whether the third set of values oscillate beyond and predetenmined threshold after a predeter- 
mined period of time and, if not. convergence has been achieved. 

19 The method of claim 1 8 wherein if no convergence has been achieved, perfbmiing the additional step of increasing 
a tim^cale for detemiining convergence beyond the predetem,ined period of time and <=07^^-f « « 
achieved, then generating a signal indicating that none of the mutations of the mutation set are present in the 

sample. 

20. The method of daim 1 5 wherein said step of determining whether a paralysis of dynamics ^^occurred oompn^s 
the steps of evaluating a Lindbald condition and, if the Lindbald condition has not been adiieved generating a 
signal indicating that paralysis of dynamics has occun^ed. 

21. The method of daim 20 wherein, if a paralysis of dynamics has occurred, performing .^^^ 
termlning whether a "mutation death" has occurred by varying a time scale for realization of the Lindbald condit«n 
and repeating the steps beginning with detemiining resonance dynamics relaxation values. 

22 The method of daim 13 wherein the step of iteratively processing the output pattern by performing a convergent 
■ reverberation also indudes the step of perfomiing a convergent reverberant dynamics ^"f"'^^^"^'!'^^.^^^ 
mutaUons using the resonance stimulus pattern to identify diagnostfc conditions represented by the mutations. 

23. The method of daim 22 wherein the step of performing a convergent revert^rent dynamics ^"^^^ 
of the mutations using the resonance stimulus pattern to identify diagnostto condrtions represented by the mutations 
indudes the steps of: 

a) determining resonance dynamics relaxation values based upon the resonance stimulus and the mutations; 

b) fi««ing the dynamics relaxation values using ensemble boundary and CSR filters to y.eld a second set of 
values* 

c) applying bulk property estimators to the dynamics relaxation values to yield a third set of values; 

d) e^luating the second and third sets of values to determine whether a predeterniined degree of resonance 
convergence has been achieved; and 

e) if no convergence, repeat steps a) - e). 

24 The method of claim 1 further including the step of filtering the diagnostic conditions identified by the convergent 
reverberant dynamics resonance analysis based upon clustering properties. 

25. The method of claim 1 wherein the biological sample is selected from a group consisting of a DNA. RNA. protein, 
peptide-nucleic acid (PNA) and targeted nucleic amplification (TNA) samples. 
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26. The method of claim 1 further including the step of rendering a diagnostic decision based on the diagnotic conditions 
identified by the convergent reverberant dynamics resonance analysis. 

27. The method of claim 26 wherein 

if said diagnostic decision is negative, then determining whether any alternatives are available; and 
if no altematives are available, selecting a new sub-set of mutations and repeating all steps. 

28. A system for analyzing an output pattern of an arrayed microdevice to identify mutations, if any. present in a 
biological sample applied to the arrayed microdevice. said system comprising: 

a tesselation unit operative to tessellate the output pattem; 

a stimulus pattem generation unit operative to generate a stimulus pattem associated with the set of known 
mutations; 

a resonance pattem generation unit operative to generate a resonance pattem representative of resonances 
between the stimulus pattem and the tessellated output pattern; 

a resonance pattem interpretation unit operative to interpret the resonance pattern to yield a set of confirmed 
mutations by comparing resonances found therein with predetemriined resonances expected for the selected 
set of mutations. 

29. The system of claim 28. further comprising: 

an extraction unit operative to extract local parametrics from the tessellated output pattem; 

a determination unit operative to detennine whether a degree of amplitude wandering representative of the 

local parametrics is within a predetermined allowable generator function limit; and 

a renormalization unit operative to renomrialize the tessellated output pattem to further match spectral prop- 
erties of the resonance pattem if the determination unit detennines that the degree of amplitude wandering is 
not within the predetermined allowable generator function limit. 

30. The system of claim 28 for performing a convergent revertserant dynamics resonance analysis of the output pattem 
to identify the mutations represented thereby, further comprising: 

a detennination unit operative to determine a first set of resonance dynamics relaxation values based upon 
the tessellated output pattem and the resonance stimulus; 

a filter operative to filter the dynamics relaxation values using ensemble boundary and CSR filters to yield a 

second set of values; . 

an estimator unit operative to apply bulk property estimators to ihe dynamics relaxation values to yield a third 

set of values; 

an evaluation unit operative to evaluate the second and third sets of values to detenmine whetiier a predeter- 
mined degree of resonance convergence has been achieved; and 

a paralysis detection unit operative to detect whether a paralysis of dynamics has occun-ed. 

31. The system according to claim 28, wherein said stimulus generation unit generates stimulus pattem based on 
Quantum Expressor Functtons. 

32. The method according to one of the claims 1 to 27 wherein an output pattem of a biochip is used as the an^yed 
microdevice output pattem. 

33. The system according to one of tiie claims 28 to 31 for analyzing an output pattem of a biochip as the arrayed 
microdevice output pattem. 



Patentanspniche 

1 . Verfahren zum Analysieren eines Ausgabemusters einer Mikroan-ayvomchtung, urn in einer auf die Mikroarrayvor- 
richtung aufgebrachten biologischen Probe, wenn ubertiaupt. vortiandene Mutationen zu identifizieren. wobei das 
Verfahren die Schritte umfasst: 

Tessellieren des Ausgabemusters. 
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Erzeugen eines der Gruppe bekannter Mutationen zugeordneten Stimutusmusters, 

Erzeugen eines Resonanzmusters, das Resonanzen zwi'schen dem Stimulusmuster und dem tessellierten 
Ausgabemuster angibt, 

Interpretieren des Resonanzmusters, um eine Gmppe bestatigter Mutationen zu eitialten, indem darin aufge- 
5 fundene Resonanzen mit fur die Gruppe von Mutationen erwarteten. vorbestimmen Resonanzen verglichen 

werden. 

2. Veifahren nach Anspruch 1 , bei dem das Ausgabemuster von einer Mikroarrayvorrichtung ertialten wird» die eine 
Mehrzahl Zellen aufweist. die jeweils eine Gruppe im Wesenttichen identischer, immobilisierter Olionukleotide auf- 

10 weisen, wobei die Olionukleotide jeder Zelle gegenuber den Olionukleotiden jeder anderen Zelle Im Wesentlichen 

eindeutig sind. 

3. Verfahren nach Anspruch 1 . bei dem der Schritt. das Ausgabemuster zu tessellieren. den Schritt umfasst, lokale 
morphologische Variationen im Ausgabemuster stark zu verstarken. 

15 

4. Verfahren nach Anspruch 3, bei dem das Stimulusmuster basierend auf Quantum-Expressor-Funktionen erzeugt 
wnrd, und bei dem der Schritt, das Ausgabemuster zu tessellieren, durchgefuhrt wird, um das Ausgabemuster so 
zu tessellieren, dass es morphologischen Charakteristika der Quantum-Expressor-Funktionen entspricht. 

20 5. Verfahren nach Anspruch 1 , femer die Schritte umfassend: 

Extrahieren lokaler Parameter aus dem tessellierten Ausgabemuster, 

Bestimmen, ob ein Ausmad einer die lokalen Parameter angebenden Amplitudenwanderung innerhalb eines 
vorbestimmten, zulassigen Generatorfunktionsgrenzwerts liegt, und 
25 wenn dies nicht der Fall ist, Renormalisieren des tessellierten Ausgabemusters, um femer mit spektralen 

Eigenschaften des Resonanzmusters ubereinzustimmen. 

6. Verfahren nach Anspaich 5, bei dem der Schritt, lokale Parameter aus dem tessellierten Ausgabemuster zu ex- 
trahieren, die Schritte umfasst, 

30 Extrahieren von Parametem, die eine integrierte Dichte von ZustSnden in jedem einer Mehrzahl von Tessel- 

lierungsbereichen in dem tessellierten Ausgabemuster angeben. 

7. Verfahren nach Anspruch 5. bei dem das AusmaB der Amplitudenwanderung bestimmt wird, indem eine Palm- 
Verteilung angewendet wird, um Generatoren zum Abschatzen einer stochastischen Wanderung zu entwickeln. 

35 

8. Verfahren nach Anspruch 1 , femer den Schritt umfassend, das tessellierte Ausgabemuster und das Stimulusmu- 
ster in ein metrisch transitives Zufallsfeld umzuwandeln. 

9. Verfahren nach Anspruch 8, femer den Schritt umfassend, das tessellierte Ausgabemuster nach einer Umwand- 
40 lung in das metrisch transitive Zufallsfeld zu renormalisieren. 

10. Verfahren nach Anspruch 9, bei dem der Renormalisierungsschritt durchgefuhrt wird, um das tessellierte Ausga- 
bemuster auf das Interval! [-n^+n] zu reskalieren. 

45 11. Verfahren nach Anspruch 1. bei dem der Schritt, ein der Gruppe bekannter Mutationen zugeordnetes Stimulus- 
muster zu erzeugen, die Schritte umfasst: 

Auswahlen einer Untergruppe der Mutationen zur Analyse, und 

Auswahlen einer Untergruppe nicht tinearer Quantum-Expressor-Funktionen aus einer Gruppe vorbestimmter, 
50 riicht linearer Quantum-Expressor-Funktionen basierend auf der ausgewdhlten Untergruppe von Mutationen. 

12. Verfahren nach Anspruch 11 , femer die Schritte umfassend: 

Transform ieren der nicht linearen Quantunrt-Expressor-Funktionen und des tessellierten Ausgabemusters in 
55 den Phasenraum. 

13. Verfahren nach Anspruch 1, bei dem der Schritt, ein Resonanzmuster zu erzeugen, den Schritt umfasst 

iteratives Verarbeiten des Ausgabemusters, indem ein konvergentes Nachhallen durchgefuhrt wird, um ein 
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Resonanzmuster zu ertialten, das Resonanzen zNMSChen einer vorgestimmten Gruppe ausgewahlter Quantum- 
Expresser-Funktionen und dem tesseilierten Ausgabemuster angibt, bis ein vorbestimmter Konvergenzgrad zwi- 
schen den in dem Resonanzmuster aufgefundenen Resonanzen und fur die Gruppe von Mutationen erwarteten 
Resonanzen erreicht wird. 

5 

14. Verfahren nach Anspruch 13, bei dem der Schritt, das Ausgabemuster iterativ zu verarbeiten, indem ein konver- 
gentes Nachhallen durchgefuhrt wird, den Schritt umfasst. eine Analyse hinsichtlich konvergenter Nachhalldyna- 
mikresonanzen des tesseilierten Ausgabemuters durchzufuhren, wobei das Resonanzstimulusmuster verwendet 
wird, um von dem tesseilierten Ausgabemuster angegebene Mutationen zu Identifizieren. 

10 

15. Verfahren nach Anspruch 14, bei dem der Schritt. eine Analyse hinsichtlich konvergenter Nachhalldynamikreso- 
nanzen durchzufuhren, die Schritte umfasst: 

a) Ermittein einer ersten Gruppe von Resonanzdynamikrelaxationswerten basierend auf dem vorverarbeiteten 
15 Ausgabemuster und dem Resonanzstimulus, 

b) Filtem der Dynamikrelaxationswerte unter Verwendung von Ensemblegrenz- und CSR-Filtem, um eine 
zweite Gruppe von Werlen zu erhalten, 

c) Anwenden von Masseeigenschaftsschatzfunktionen auf die Dynamikrelaxationswerte, um eine dritte Grup- 
pe von Werten zu erhalten, 

20 d) Evaluieren der zweiten und dritten Gruppe von Werten, um einen Resonanzkonvergenzgrad zu ermittein. 

und 

e) Ennittein aus dem Resonanzkonvergenzgrad. ob eine Dynamiksperrung aufgetreten ist. und, wenn dies 
der Fall ist. wiederholen der Schritte a) - e). 

25 16. Verfahren nach Anspruch 15, bei dem der Schritt, basierend auf dem vorverarbeiteten Ausgabemuster und dem 
Resonanzstimulus Resonanzdynamikrelaxationswerte zu bestimmen, die Schritte umfasst 

Anwenden eines faltungslosen Entwicklungsoperators auf das tessellierte Ausgabemuster. 

17. Verfahren nach Anspruch 15, bei dem der Schritt, Masseeigenschaftsschatzfunktionen auf die Dynamikrelaxati- 
30 onswerte anzuwenden, um eine dritte Gruppe von Werten zu erhalten, die Schritte umfasst 

Anwenden eines Kopplungsoperators, mit dem die das tessellierte Ausgabemuster angebenden Dynamikre- 
laxationswerte gekoppelt sind, auf einen nicht linearen Informationsfilter. 

18. Verfahren nach Anspruch 15. bei dem der Schritt, die zweiten und dritten Gruppen von Werten zu evaluieren, um 
35 einen Resonanzkonvergenzgrad zu ermittein. die Schritte umfasst 

Ermittein, ob die dritte Gruppe von Werten nach einer vorbestimmten Zeitdauer uber einen vorbestimmten 
Grenzwert hrnaus oszillieren und, falls dies nicht der Fall ist, Konvergenz en-eicht wurde. 

19. Verfahren nach Anspmch 18, bei dem, wenn keine Konvergenz erreicht wurde, derzusatzliche Schritt durchgefuhrt 
40 wird, einen ZeltmaRstab zu vergroBem, um nach der vorbestimmten Zeitdauer Konvergenz festzustellen, und, 

wenn Konvergenz weiterhin nicht en-eicht wird, dann ein Signal zu erzeugen. das angibt. dass keine Mutationen 
der Mutationsgruppe in der Probe voriiegen. 

20. Verfahren nach Anspruch 1 5, bei dem der Schritt, zu emnittein, ob eine Dynamiksperrung aufgetreten ist, die Schrit- 
45 te umfasst, eine Lindbald-Bedingung zu evaluieren, und, wenn die Lindbald-Bedingung nicht erreicht wurde. ein 

Signal zu erzeugen. das angibt, dass eine Dynamiksperrung aufgetreten ist. 

21 . Verfahren nach Anspruch 20. bei dem, wenn eine Dynamiksperrung aufgetreten ist, die zusatzlichen Schritte aus- 
gefOhrt werden. festzustellen. ob ein "Mutationstod" aufgetreten ist. indem ein ZeitmaRstab zur Realisierung der 

50 Lindbald-Bedingung variiert und die Schritte beginnend mit der Ermittlung von Dynamikresonanzrelaxationswerten 

wiederholt werden. 

22. Verfahren nach Anspmch 1 3. bei dem der Schritt. das Ausgabemuster iterativ zu verart^eiten, indem ein konver- 
gentes Nachhallen durchgefuhrt wird, auch den Schritt umfasst, eine Analyse hinsichtlich konvergenter Nachhall- 

55 dynamikresonanzen der Mutationen durchzufuhren, wobei das Resonanzstimulusmuster verwendet wird, um von 

den Mutationen angegebene. diagnostische Zustande zu identifizieren. 

23. Verfahren nach Anspruch 22. bei dem der Schritt, eine Analyse hinsichtlich konvergenter Nachhalldynamikreso- 
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nanzen der Mutationen unter Verwendung des Resonanzstimutusmusters durchzufOhren, um von den Mutationen 
angegebene, diagnostische ZustSnde zu identifizieren, die Schritte umfasst: 

a) Ermitteln einer ersten Gmppe von Resonanzdynamikrelaxationswerten basierend auf dem Resonanzsti- 
5 mulus und den Mutationen, 

b) Filtem der Dynamlkrelaxationswerte unter Verwendung von Ensemblegrenz- und CSR-Filtem. um eine 
zweite Gruppe von Werten zu erhalten, 

c) Anwenden von MasseeigenschaftsschStzfunktionen auf die Dynamikrelaxationswerte, um eine dritte Grup- 
pe von Werten zu ertialten, 

10 d) Evaluieren der zwetten und dritten Gruppe von Werten, um einen Resonanzkonvergenzgrad zu enmittein, 

und 

e) Ermitteln aus dem Resonanzkonvergenzgrad. ob eine Dynamiksperrung aufgetreten ist, und, wenn dies 
der Fall ist. wlederholen der Schritte a) - e). 

15 24. Verfahren nach Anspruch 1 , femer den Schritt umfassend, die durch die Analyse hinsichtlich konvergenter Nactv 
halldynamikresonanzen identifizierten diagnostischen Zustande basierend auf Clustereigenschaften zu filtem. 

25. Verfahren nach Anspruch 1, bei dem die biologische Probe aus einer Gruppe ausgewahit ist, die aus DNA-, RNA-, 
Protein-, Peptid-Nucleinsaure-(PNA)- und ZielzelIkemverstarkungs-(TNA)-Proben besteht. 

20 

26. Verfahren nach Anspruch 1, femer den Schritt umfassend. basierend auf den durch die Analyse hinsichtlich kon- 
vergenter Nachhalldynamikresonanzen identifizierten, diagnostischen ZustSnden eine diagnostische Entschei- 
dung zu treffen. 

25 27. Verfahren nach Anspruch 26, bei dem 

wenn die diagnostische Entscheidung negativ ist, dann ermittelt wird, ob irgendwelche Altemativen verfugbar 

sind, und 

wenn keine Alternativen verfugbar sind, eine neue Untergruppe von Mutationen ausgewdhit und al!e Schritte 
wiederholt werden. 

30 

28. System zum Analysieren eines Ausgabemusters einer Mikroarrayvorrichtung. um in einer auf die Mikroarrayvor- 
richtung aufgebrachten biologischen Probe, wenn Qt>erhaupt, vorhandene Mutationen zu identifizieren, wobei das 
System umfasst: 

35 eine Tessellierungseinheit. die betriebsfShig ist. das Ausgabemuster zu tessellieren, 

eine Stimulusmustererzeugungseinheit. die betriebsfdhig ist, ein Stimulusmuster zu erzeugen, das der Gruppe 
bekannter Mutationen zugeordnet ist, 

eine Resonanzmustererzeugungseinheit, die betriebsfahig ist, ein Resonanzmuster zu erzeugen, das Reso- 
nanzen zwischen dem Stimulusmuster und dem tessellierten Ausgabemuster angibt, 
40 eine Resonanzmusterlnterpretationseinheit, die betriebsfahig ist, das Resonanzmuster zu interpretieren, um 

eine Gruppe bestatigter Mutationen zu erhalten, indem darin aufgefundene Resonanzen mit fur die Gmppe 
von Mutationen erwarteten, vorbestimmten Resonanzen verglichen werden. 

29, System nach Anspnjch 28, femer umfassend: 

45 

eine Extraktionseinheit, die betriebsfahig ist, lokale Parameter aus dem tessellierten Ausgabemuster zu ex- 
trahieren, 

eine Ermittlungseinheit, die betriebsfdhig ist, zu ermitteln, ob ein AusmaH einer die lokalen Parameter ange- 
benden Amplitudenwanderung in einem vorbestimmten, zulSssigen Generatorfunktlonsgrenzwert liegt, und 
50 eine Renormalisierungseinheit, die betriebsfdhig ist, das tessellierte Ausgabemuster zu renormalisieren, um 

femer mit spektralen Eigenschaften des Resonanzmusters ubereinzustimmen, wenn die Ermittlungseinheit 
feststellt, dass das Ausmali der Amplitudenwandemng nicht innerhalb des vorgestimmten, zulassigen Gene- 
ratorfunktionsgrenzwerts liegt. 

55 30. System nach Anspruch 28, um eine Analyse hinsichtlich konvergenter Nachhalldynamikresonanzen des Ausga- 
bemusters durchzufOhren, um von diesem angegebene Mutationen zu identifizieren, femer umfassend: 

eine Ermittlungseinheit. die betriebsfahig ist, eine erste Gruppe von Resonanzdynamikrefaxationswerten ba- 
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sierend auf dem tessellierten Ausgabemuster und dem Resonanzstimulus zu bestimmen, 

einen Filter, der betriebsfahig ist, die Dynamil<relaxationswerte unter Verwendung von Ensemblegrenz- und 

CSR-Filtem zu filtem, um eine zweite Gruppe von Werten zu erhalten. 

eine SchSlzeinheit, die betriebsfahig ist, Masseeigenschaflsschatzfunktionenauf die Dynamikrelaxationswer- 
5 te anzuwenden, um eine dritte Gruppe von Werten zu ertialten, 

eine Evaluiemngselnheit. die betriebsfahig ist, die zweiten und dritten Gmppen von Werten zu evaluieren, um 
zu ermitteln, ob ein vorbestimmter Resonanzkonvergenzgrad erreicht wurde. und 

eine Spenrungsdetektionseinheit, die betriebsfahig ist, zu ermitteln. ob eine Dynamiksperrung aufgetreten ist. 

10 31 . System nach Anspmch 28. bei dem die Stimuluserzeugungseinheit basierend auf Quantum-Expressor-Funktionen 
ein Stimulusmuster erzeugt. 

32. Verfahren nach einem der Anspruche 1 bis 27, bei dem ein Ausgabemuster eines Biochips als das Ausgabemuster 
einer Mikroarrayvorrichtung verwendet wird. 

15 

33. System nach einem der Anspruche 28 bis 31 zum Anaiysieren eines Ausgabemusters eines Biochips als das 
Ausgabemuster einer Mikroarrayvorrichtung. 

20 Revendications 

1. Procede d*analyse d'un motif de sortie d'un mrcrodispositif ordonn6 pour identifier des mutations, le cas ech^ant, 
pr^sentes dans un echantillon biologique applique au microdispositif ordonn^. ledit proc6d§ comprenant les stapes 
consistant d : 

25 

paver le motif de sortie ; 

generer un motif de stimulus associe a un ensemble de mutations connues ; 

gen6rer un motif de resonance representatrf de resonances entre le motif de stimulus et le motif de sortie pave ; 
interpreter le motif de resonance pour produire un ensemble de mutations confirmees en comparant les r6- 
30 sonances decouveries dans celui-ci a des resonances pr6determinees attendues pour I'ensemble setectionnS 

de mutations. 

2. Procede de la revendication 1 . dans lequel le motif de sortie est issu d*un microdispositif ordonne ayant une pluralite 
de cellules ayant chacune un ensemble d'oligonudeotides quasiment identiques immobilises avec les oligonu- 

35 deotides de chaque cellule etant quasiment uniques par rapport aux oligonudeotides de chaque autre cellule. 

3. Procede de la revendication 1 , dans lequel ladite etape de pavage du motif de sortie comprend I'etape consistant 
a amplifier des variations morphologiques fortement locales dans le motif de sortie. 

40 4. Procede de la revendication 3, dans lequel le motif de stimulus est genere a partir de fonctions d'expression 
quantique et dans lequel ladite etape de pavage du motif de sortie est executee pour paver le motif de sortie afin 
qu'il corresponde d des caracteristiques morphologiques des fonctions d'expression quantique. 
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5. Procede de la revendication 1 comportant en outre les etapes consistant d : 



extraire une parametrique locale du motif de sortie pave ; 

determiner si un degre de variation d'amplitude representatif de la parametrique locale se situe e Tinterieur 
d'une limite de fonction generatrice autorisee predeterminee ; et 

sinon. renormaliser le motif de sortie pave pour qu'il con-esponde davantage e des proprietes spectrales du 
so motif de resonance. 

6. Precede de la revendication 5, dans lequel ladite etape d'extraction d'une parametrique locale du motif de sortie 
pave comprend retape consistant a extraire des param6tres representatifs d'une densite d'etats integree au sein 
de chacune d'une pluralite de regions de pavage d I'interieur du motif de sortie pave. 



Procede de la revendication 5, dans lequel un degre de variation d'amplitude est determine en appliquant une 
distribution de Palm afin de developper des generateurs pour estimer une variation stochastique. 
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8. Proc6d6 de la revendication 1 oomportant en outre l*§tape consistant d transformer le motif de sortie pave et le 
motif de stimulus en un champ al^atoire m^triquement transitif. 

9. Proc6d6 de la revendication 8 comportant en outre I'^tape consistant d renormaliser le motif de sortie pav6 aprds 
5 transfonmation en champ al6atoire m6triquement transitif. 

10. Proc6d§ de la revendication 9. dans lequel ladite 6tape de renomialisation est execut^e pour r^duire le motif de 
sortie pav6 d Tintervalle [-re, +rel. 

10 11. Proems de la revendication 1, dans lequel I'Stape de production d'un motif de stimulus associS a un ensemble 
de mutations connues comprend les etapes consistant d : 

s6lectionner un sous-ensemble de mutations pour analyse ; et 

s6lectionner un sous-ensemble de fonctions d'expression quantique non lineaires a partir d'un ensemble de 
^5 fonctions d'expression quantique non lineaires predetemnrnees sur la base du sous-ensemble de mutations 

selectionn6. 

12. Proced6 de la revendication 11 comportant en outre I'etape consistant ^ transformer les fonctions d'expression 
quantique non lineaires et le motif de sortie pavS en espace des phases. 

20 

13. Proced§ de la revendication 1, dans lequel I'etape de production d'un motif de resonance comporte T^tape con- 
sistant a traiter par iteration le motif de sortie en executant une reverberation convergente afin de produlre un motif 
de resonance repr^sentatif de resonances entre un ensemble predetermine de fonctions d'expression quantique 
seiectionnees et le motif de sortie pave jusqu'e ce qu'un degre predetermine de convergence soit attaint entre les 

25 resonances decouvertes dans le motif de resonance et des resonances attendues pour I'ensemble de mutations. 

14. Procede de la revendication 13, dans lequel I'etape de traltement iteratif du motif de sortie par execution d'une 
reverberation convergente comporte I'etape consistant a executer une analyse de resonance par dynamlque re- 
verber6e convergente du motif de sortie pave a I'aide du motif de stimulus de resonance afin d'identifier des 

30 mutations representees par le motif de sortie pave. 

15. Precede de la revendication 14, dans lequel I'etape d'execution d'une analyse de resonance par dynamique re- 
verberee convergente comporte les etapes consistant d : 

35 a) determiner un premier ensemble de valours de relaxation de la dynamique de resonance a partir du motif 

de sortie preconditionne et du stimulus de resonance ; 

b) filtrer les valeurs de relaxation de la dynamique d I'aide de filtres de limite d'ensemble et CSR afin de 
produire un second ensemble de valeurs ; 

c) appliquer des estimateurs de proprietes grossiers aux valeurs de relaxation de la dynamique pour produire 
^0 un troisieme ensemble de valeurs ; 

d) evaluer les second et troisieme ensembles de valeurs pour determiner un degre de convergence de 
resonance ; et 

e) determiner e partir du degre de convergence de resonance si une paralysie de la dynamique s'est produite 
et. si c'est le cas, repeter les etapes a) - e). 

45 

1 6. Procede de la revendication 1 5, dans lequel ladite etape de determination de valeurs de relaxation de la dynamique 
de resonance ^ partir du motif de sortie preconditionne et du stimulus de resonance comprend I'etape consistant 
d appliquer un operateur d'evolution sans convolution au motif de sortie pave. 

50 17. Procede de la revendication 15, dans lequel ladite etape d'application d'estimateurs de proprietes grossiers aux 
valeurs de relaxation de la dynamique pour produire un troisieme ensemble de valeurs comprend I'etape consistant 
d appliquer un operateur de couplage qui couple les valeurs de relaxation de la dynamique representatives du 
motif de sortie pave e un filtre d'information non lineaire. 

55 18. Precede de la revendication 15, dans lequel ladite etape d'evaluatlon des second et troisieme ensembles de va- 
leurs pour determiner un degre de convergence de resonance comprend I'etape consistant a determiner si le 
troisieme ensemble de valeurs oscille au-del^ d'un seuil predetermine apres une period e de temps predetenninee 
et, sinon. si une convergence a ete atteinte. 
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19. Proc6d6 de la revendication 18. dans lequel si aucune convergence n'a 6t§ altefnte, on execute T^tape supp!6- 
mentaire consistant a augmenter r^cheile de temps pour determiner une convergence au-deld de la p^riode de 
temps pr^d^terminee et, si une convergence n'est toujours pas atteinte, on g^nere alors un signal indiquant 
qu*aucune des mutations de Tensemble de mutations n*est pr^sente dans I'^chantillon. 

5 

20. Proc^d^ de la revendication 15, dans lequel ladite dtape consistant a d^temniner si une paralysie de la dynamique 
s'est produite comprend les stapes consistant ^ ^valuer une condition de Lindbald et. si la condition de Lindbald 
n*a pas 6t6 atteinte. genSrer un signal indiquant qu*une paralysie de la dynamique s'est produite. 

10 21 . Proc6d6 de la revendication 20, dans lequel, si une paralysie de la dynamique s'est produite, on execute les Stapes 
supplemental res consistant d determiner si une "mort des mutations" s'est produite en faisant varier rechelle de 
temps pour la realisation de la condition de Lindbald et on repete les etapes en commengant par determiner les 
valeurs de relaxation de la dynamique de resonance. 

15 22. Precede de la revendication 13, dans lequel I'etape de traitement Iteratif du motif de sortie par execution d'une 
reverberation convergente comporte egaiement I'etape consistant a executer une analyse de resonance par dy- 
namique reverberee convergente des mutations a I'aide du motif de stimulus de resonance afin d'identifier des 
conditions de diagnostic representees par les mutations. 

20 23. Precede de la revendication 22. dans lequel i'etape d'execution d'une analyse de resonance par dynamique re- 
verberee convergente des mutations a I'aide du motif de stimulus de resonance afin d'identifier des conditions de 
diagnostic representees par les mutations comporte les etapes consistant d : 

a) determiner un premier ensemble de valeurs de relaxation de la dynamique de resonance a partirdu stimulus 
25 de resonance et des mutations ; 

b) fittrer les valeurs de relaxation de la dynamique d I'aide de filtres de limite d'ensemble et CSR pour produire 
un second ensemble de valeurs ; 

c) appliquer des estimateurs de proprietes grossters aux valeurs de relaxation de la dynamique pour produire 

un troisieme ensemble de valeurs ; 
30 d) evaluer les second et troisieme ensembles de valeurs pour determiner si un degre predetermine de con- 

vergence de resonance a ete atteint ; et e) en cas d'absence de convergence, repeter les etapes a) - e). 

24. Precede de la revendication 1 comportant en outre i'etape consistant a filtrer les conditions de diagnostic identifiees 
par I'analyse de resonance par dynamique reverberee convergente d partir de proprietes de groupage. 
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25. Precede de la revendication 1 . dans lequel rechantlllon biologique est choisi dans un groupe constitue par des 
echantillons d'ADN, d'ARN. de proteine. d'acide peptide-nudeique (PNA) et d'amplification nucieique cibiee (TNA). 

26. Precede de la revendication 1 comportant en outre i'etape consistant a rendre une decision de diagnostic d partir 
40 des conditions de diagnostic identifiees par I'analyse de resonance par dynamique reverberee convergente. 

27. Procede de la revendication 26. dans lequel. 

si ladite decision de diagnostic est negative, on determine alors si une quelconque alternative est disponibte ; 

et 

^5 si aucune alternative n'est disponibie. on seiectionne un nouveau sous-ensemble de mutations et on repete 

toutes les etapes. 

28. Systeme d'analyse d'un motif de sortie d'un microdispositif ordonne pour identifier des mutations, le cas echeant. 
presentes dans un echantillon biologique applique au microdispositif ordonne, ledit systeme comprenant : 

50 

une unite de pavage operant pour paver le motif de sortie ; 

une unite de production de motif de stimulus operant pour generer un motif de stimulus associe d I'ensembie 
de mutations connues ; 

une unite de production de motif de resonance operant pour generer un motif de resonance representatif de 
55 resonances entre un motif de stimulus et le motif de sortie pav6 ; 

une unite d'interpretatlon de motif de resonance operant pour interpreter le motif de resonance afin de produire 
un ensemble de mutations confirmees en comparant des resonances decouvertes dans cetui-ci e des reso- 
nances predetenminees attendues pour I'ensembie seiectionne de mutations. 
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29. Syst6me de la revendication 28, comprenant en outre : 

une unit6 d'extraction operant pour extraire une param6trique locale du motif de sortie pav6 ; 
une unit6 de determination operant pour determiner si un degr6 de variation d'amplitude repr6sentatif de la 
5 parametrique locale se situe d I'interieur d'une limite de fonction g6n6ratrice autoris^e pr6d6termin6e ; et 

une unite de renormalisation operant pour renonmaliser le motif de sortie pave afin qu'il con-esponde davantage 
d des proprietes spectrales du motif de resonance si Tunite de detemnination determine que le degre de va- 
riation d'amplitude ne se situe pas d I'interieur de la limite de fonction generatrice autorisee predeterminee. 

10 30. Systeme de la revendication 28 pour ex6cuter une analyse de resonance par dynamique r6vert5er6e convergente 
du motif de sortie afin tf identifier les mutations representees par ce!ui-ci. comprenant en outre : 

une unite de determination operant pour determiner un premier ensemble de valeurs de relaxation de ia dy- 
namique de resonance ^ partir du motif de sortie pave et du stimulus de resonance ; 
IS un filtre operant pour filtrer les valeurs de relaxation de la dynamique d I'aide de filtres de limite d'ensemble 

et CSR afin de produire un second ensemble de valeurs ; 

une unite estimatrice operant pour appliquer des estimateurs de proprietes grossiers aux valeurs de relaxation 
de la dynamique afin de produire un troisi6me ensemble de valeurs ; 

une unite d'6valuation operant pour evaluer les second et troisieme ensembles de valeurs afin de detemiiner 
20 si un degre predetermine de convergence de resonance a ete attaint ; et 

une unite de detection de paralysie operant pour detecter si une paralysie de la dynamique s'est produite. 

31. Systeme selon la revendication 28. dans lequel ladite unite de production de stimulus genere un motif de stimulus 
d partir de fonctions d'expression quantique. 

25 

32. Precede selon Tune des revendications 1 a 27, dans lequel un motif de sortie tfune puce S ADN est utilise comme 
le motif de sortie du microdispositif ordonne. 

33. Systeme selon I'une des revendications 28 d 31 pour analyser un motif de sortie d'une puce d ADN comme le 
30 motif de sortie du microdispositif ordonne. 
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